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Разработана математическая модель ректификации формалина-сырца, учитывающая химические превращения и позволяющая 
проводить численный анализ эффективности разделения компонентов технологических потоков при различных режимах рабо- 
ты массообменной колонны и при изменении состава формалина-сырца. Выявлено различие в схемах взаимодействия фор- 
мальдегида с водой и метанолом. Предложен механизм пересчета аналитических концентраций компонентов и расходов техно- 
логических потоков в составы и расходы псевдобинарных потоков. Установлены интервалы изменения межтарельчатого уноса 
при различных режимах орошения колонны и показана необходимость принятия мер по увеличению эффективности действую- 
щих контактных устройств. 


Формалин, представляющий собой водный ра- 
створ формальдегида СН 2 0 и метанола СН 3 ОН, на 
сегодняшний день является одним из важнейших 
продуктов нефтехимической промышленности, 
спрос на который продолжает расти. По причине 
пожаро- и взрывоопасности предприятий данной 
отрасли методы математического моделирования 
более привлекательны по сравнению с методами 
моделирования физического, если требуется прог- 
ноз производительности промышленной установ- 
ки или поиск оптимального режима ее функциони- 
рования. Математическая модель ректификации 
формалина-сырца позволяет получить профиль 
концентраций разделяемых веществ по высоте ко- 
лонны и выполнить количественный анализ эф- 
фективности разделения при различных техноло- 
гических режимах и составах формалина-сырца. 


На сегодняшний день существует ряд подходов 
к моделированию процесса в ректификационных 
колоннах. В основе структуры любой модели рек- 
тификации лежат уравнения, описывающие меж- 
фазный перенос компонентов разделяемой смеси. 
Эти уравнения позволяют определить состав паро- 
вой фазы, покидающей ступень разделения, по из- 
вестному составу жидкости на ней и составу посту- 
пающего пара. Математическое описание межфаз- 
ного равновесия осложняется, в основном, много- 
компонентностью и неидеальностью растворов. 
Водно-метанольный раствор формальдегида, обра- 
зующийся в результате взаимодействия формаль- 
дегида, метанола и воды (система Ф-М-В), отно- 
сится к многокомпонентным, неидеальным систе- 
мам, осложненным химическими взаимодействи- 
ями [ 1 ] . Традиционно точность модели, в том числе 
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и модели парожидкостного равновесия, устанавли- 
вается путем сравнения результатов расчета с экс- 
периментальными данными. Как правило, экспе- 
риментально определяют концентрации индивиду- 
альных компонентов, образующих раствор, в жид- 
кой и паровой фазах. При таком определении мо- 
жет быть утрачена информация об истинном ком- 
понентном составе раствора, обусловленном хими- 
ческими реакциями между образующими его веще- 
ствами и, следовательно, о коэффициентах актив- 
ности соединений, существующих только в раство- 
ре. Система Ф-М-В не является исключением, т.е. 
ее равновесные характеристики (равновесные со- 
ставы фаз и коэффициенты активности), получен- 
ные экспериментальным путем, никак не связаны 
с качественным составом раствора. Существующие 
на сегодняшний день математические модели па- 
рожидкостного равновесия этого раствора ориен- 
тированы на оценку коэффициентов активности 
только формальдегида, метанола и воды. 

Учитывая способность формальдегида реагиро- 
вать с водой и метанолом параллельно и независимо 
в смешанном растворителе [2] нами предложена со- 
вокупная схема [3], основанная на реакциях, которые 
протекают в бинарных растворах (Ф-М, Ф-В) [1,4]. 

Ранее в основе математического описания рав- 
новесных свойств данной системы как содержащей 
реакции полимеризации использовалась схема, ос- 
нованная на представлении об аналогии реакций 
формальдегида с водой и реакций формальдегида с 
метанолом [5]. Как видно из табл. 1, такая аналогия 
отсутствует. 


Таблица 1 . Схемы реакций формальдегида в бинарных ра- 
створах 


Водный раствор формальдегида 

Метанольный раствор 
формальдегида 

Р+А<-»РА (1) 

Р+Мо-РМ (II) 

РА+РА<->Р 2 А+А 

Р+РМ<->Р 2 М 

Р 2 А+РА<-»Р 2 А+А 

Р+Р 2 М<->Р 3 М 

Р гг1 А+ РА<-»Р П А+ А 

Р+Рпг,М<->Р т М 


Здесь С ~ формальдегид; М ~ метанол; А ~ вода; Р,А (і=[1,п]) 
- полиоксиметиленгидраты; Р/М (]=[1,т]) - полиоксимети- 
ленсольваты; п, гл - максимально возможное количество ато- 
мов формальдегида в молекуле полиоксиметиленгидрата и 
полиоксиметиленсольвата, соответственно 

Теория ассоциированных растворов и математи- 
ческое описание фазового равновесия систем, ком- 
поненты которых вступают в химические взаимо- 
действия - это основные инструменты установле- 
ния взаимосвязи между коэффициентами активно- 
сти компонентов, из которых образован раствор, и 
содержанием молекул этих веществ в растворе. Ис- 
пользование перечисленных инструментов приме- 
нительно к системе Ф-М-В сопряжено с необходи- 
мостью, во-первых, точной оценки качественного и 
количественного состава раствора и, во-вторых, с 
выбором наиболее значимых превращений. 

Анализ многочисленных исследований компо- 
нентного состава В-М растворов формальдегида 
показал, что в промышленной ректификационной 


колонне разделяемая смесь наряду с метанолом и 
водой может содержать более 10 веществ [6]. Воз- 
можность установления конкретной совокупности 
реакций в растворах формальдегида, и, как след- 
ствие, компонентного состава раствора при тех или 
иных условиях и на сегодняшний день представля- 
ет собой самостоятельную научную задачу. 

Важность учета влияния химических реакций под- 
тверждают существенные различия между коэффици- 
ентами активности, определенными из эксперимен- 
тальных данных, и рассчитанными одним из наибо- 
лее универсальных методов, с учетом влияния струк- 
туры только формальдегида, метанола и воды [4]. 

Как известно, при ректификации происходит 
разделение раствора на индивидуальные компо- 
ненты или их смеси, различающиеся температура- 
ми кипения. О наличии продуктов химических ре- 
акций в системе Ф-М-В свидетельствует тот факт, 
что, несмотря на то, что минимальной температу- 
рой кипения обладает формальдегид (его темпера- 
тура кипения -19 °С), в колонне в качестве дистил- 
лята отбирается почти чистый метанол (его темпе- 
ратура кипения 64,5 °С). 

Кроме того, в дистилляте отсутствует формальде- 
гид, что позволяет нам утверждать, что он полностью 
расходуется на реакции с метанолом, водой и по- 
лиоксиметилсольватами. Нижние пределы интерва- 
лов экспериментально полученных значений темпе- 
ратур кипения растворов формальдегида (у метаноль- 
ного раствора [66, 6.. .80,0]; у водного [99,1. ..100,0]; у 
водно-метанольного [70,1. ..95, 8] при атмосферном 
давлении [7, 8]) существенно выше температуры ки- 
пения формальдегида, что также свидетельствует о 
его отсутствии в растворе, поскольку теоретический 
интервал температур кипения водных растворов дол- 
жен находиться в пределах [-19,0.. .64,5], а метаноль- 
ных - в интервале [-19,0.. .100,0]. 

Поскольку в ряде работ [1-3, 7, 8], посвящен- 
ных экспериментальному изучению парожидко- 
стного равновесия в данной системе, допущение об 
отсутствии паров формальдегида и продуктов реак- 
ции в паровой фазе признано уместным, коэффи- 
циент активности низкокипящего компонента 
(НК) в разработанной модели процесса рассчиты- 
вается по следующему уравнению: 

Р 

* “ к р° . ѵ +р° П-Х Г 

здесь Р— общее давление в системе, Па; Р° - давле- 
ние паров чистого вещества, Па; х - мольная доля 
вещества в жидкой фазе. 

Для моделирования ректификации формали- 
на-сырца система Ф-М-В, находящаяся в колонне, 
рассматривается как псевдобинарный раствор, под 
низкокипящим компонентом которого подразуме- 
вается метанол. Поскольку присутствие продуктов 
реакций наряду с водой в составе раствора бес- 
спорно, второй, высококипящий, компонент (ВК) 
представляет собой остальную часть раствора. На 
основе значений констант равновесия [4] реак- 
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ций I, II (табл. 1), в которые вступает формальдегид 
в смешанном растворителе, расходы и составы по- 
токов в ректификационной колонне с учетом про- 
дуктов химических реакций рассчитываются по ал- 
горитму, приведенному ниже. 


здесь С - расход технологического потока, кг/ч. 

3. Расчет расходов псевдобинарных потоков С, 
кмоль/ч: 

С'=Х(С 7 /Мг.), 


Алгоритм расчета составов 

1. Определение числа молей метиленгидратов 
п м выполняется на основе соотношения констант 
равновесия [2, 8] первых реакций (1, 2), в которые 
вступает формальдегид в смешанном растворителе: 

СН 2 0+Н 2 0 < ''" ' - х СНЛОН). (1) 

СН 2 0+СН,0Н < А >гм _ > сНзОСН ,ОН (2) 

„ _ Х Р ' Е РрЛ 

П РА ТГ , Т г ’ 

РРрА + К РрМ 

здесь х г - аналитическая концентрация формаль- 
дегида (мол. %). 

2. Определение числа молей метиленсольватов п ы : 

х р ' КРрм 

П РМ тг , т , 

РР РА + КРрм 

3. Расчет числа молей воды, не вступившей в 
реакции А г '“, как остатка от исходного (аналитиче- 
ского) ее содержания на 100 моль раствора: 

П л #ее — Х А~ П РА ' 

4. Расчет числа молей метанола, не вступивше- 
го в реакции М ге ‘, как остатка от исходного (анали- 
тического) его содержания на 100 моль раствора: 

П м Р“ = Х м — П РМ • 

5. Определение мольных долей х ГА , х ш , х/™, х^\ 

х, =100- и ( / ^ и, , 

і 

где / - РА, РМ, Л г “, Ш“. 

6. Определение мольных долей А / ' ее (х / ’р™) и 
М'“(х^гее) как компонентов псевдобинарной системы: 

- суммирование х,/™ и х/«: 

- вычисление: 

х? = 100-х,./ У/ , 

здесь х у - аналитическая концентрация метанола 
или воды, соответственно; ]-Л гее ,М г ‘ е . 

Алгоритм расчета расходов (кмоль/ч) 
псевдобинарных потоков 

1. Расчет массовых долей РА, РМ, АР”, М'“\ 

Мг -х. 

Х ) — \ д , ’ 

2 _а Мг і- х і 

і 

здесь РА, РМ, АР“, М ке , Мг, - молекулярная масса 
вещества / 

2. Расчет расходов потоков «по веществу» О р кг/ч: 

С ] =~х ] -С, 


здесь і=А ке , М гее . 

Результаты пересчета составов и расходов реаль- 
ных потоков в псевдобинарные с учетом химических 
реакций в потоках колонны представлены в табл. 2. 


Таблица 2. Результаты пересчета составов и расходов ректи- 
фикационной колонны 


Наименова- 
ние потока 

] Содержание веществ \ 

\ мол. % 

мае. % | 

РА 

РМ 


М*' 

РА 

РМ 

Л” 

М*" 

Метанол- 

ректификат 

0,41 

0,05 

3,71 

95,83 

0,62 

0,11 

2,12 

97,16 

Формалин- 

сырец 

25,73 

3,40 

48,43 

22,44 

40,68 

6,95 

28,72 

23,65 

Формалин 

товарный 

32,72 

4,33 

60,72 

2,24 

52,29 

8,93 

36,39 

2,38 

Параметры псевдобинарных потоков | 

Наименова- 
ние потока 

! расход 

состав, мол. % | 

кг/с 

кмоль/ч 

АР* 

М"' 

Метанол- 

ректификат 

1,20 

0,04 

3,73 

96,27 

Формалин- 

сырец 

2,82 

0,07 

68,34 

31,66 

Формалин 

товарный 

1,62 

0,03 

96,45 

3,55 


Расчеты по разработанной модели, основанной 
на методе счета «от тарелки к тарелке», выполня- 
ются в направлении от куба к дефлегматору по сле- 
дующему алгоритму. 

1. Расчет расходов и составов псевдобинарных 
технологических потоков. 

2. Расчет изменения давления (Р, Па) и темпера- 
туры ( Т, К) по высоте колонны как функции да- 
вления по аппроксимирующей зависимости: 

Т = 595,001-е _4О522 ' 39/р . 


3. Расчет давления паров чистых веществ (воды и 
метанола) в зависимости от температуры и да- 
вления по высоте колонны методом Риделя- 
Планка-Миллера. 

4. Исчерпывающая часть ректификационной ко- 
лонны: 


4.1. 


Расчет составов жидкости х, и пара у, на 
первой тарелке: 

_ (! 'У:< , • " ; 


У 1,1 =У № .і+ е - 


IV, і 




і \ѵ ,у 


•*.,/ -Уи 


« 1 


здесь ] - компонент технологического потока; С, 
]Ѵ- расход потока пара и кубового остатка, соответ- 
ственно, кмоль/ч; Р, Ф - расход питания и ороше- 
ния колонны, соответственно, кмоль/ч, Хщ - содер- 
жание в жидкой фазе кубового остатка компонента у, 
мольные доли; у ѴІ - содержание в паровой фазе ку- 
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бового остатка компонента у, мольные доли; - ко- 

эффициент активности компонента у в кубовом ос- 
татке; Е— эффективность контактного устройства. 

4.2. Расчет составов жидкости и пара в исчер- 
пывающей части колонны, начиная со 2 
тарелки: 


\і=- 


С-У,_ Х :+Ш-Х , 




Р + Ф 


У и =Уі 


І-1.У 


+ Е- 


Р. 

Уі.і -.'м../ 


здесь і - порядковый номер контактного устройства. 

5. Укрепляющая часть ректификационной колонны. 
5.1. Расчет составов жидкости и пара в укре- 
пляющей части: 


‘І’У./ • »*.,./ • ('У ,./ 


X, , - 



Ф + Р 


( ^ 

= Уі-щ + е ■ 

Уі.і ^- Х і.і-Ушт 


V У 


здесь Хц - содержание в дистилляте компонента у, 
мольные доли; Р, - давление на і- ой тарелке. 

6. Расчет массовых концентраций НК и ВК на ос- 
нове мольных, полученных по пп. 4, 5. 

7. По высоте колонны рассчитываются меняю- 
щиеся с изменением состава физические свой- 
ства потоков: 

7.1. плотность парового потока ; 

7.2. плотность потока жидкости ; 

7.3. поверхностное натяжение потока жидкости. 

8. Расчет эффективности колпачковой тарелки по 
высоте колонны. 

Допущения в модели: 

1 . Истинный компонентный состав системы Ф-М-В 
в колонне: Н 2 0; СН 2 (ОН) 2 ; СН 3 ОН; СН 3 ОСН 2 ОН. 

2. Исходная смесь рассматривается как псевдоби- 
нарная, низкокипящим компонентом которой 
является метанол. 

3 . Питание подается при температуре кипения смеси. 

4. Количество жидкости (нагрузка по жидкости) 
по высоте секции постоянна. 

5. Количество поднимающегося по колонне пара 
постоянно по высоте колонны. 

6. Флегма на орошение подается при температуре 
кипения. 

7. Колонна снабжена полным дефлегматором. 

8. Кипятильник работает как полный испаритель. 

9. Массопередача на тарелках колонны эквимо- 
лярная. 

10. В зоне массобмена на тарелке осуществляется 
полное вытеснение по пару. 

Проверка целесообразности разработанной мо- 
дели выполнена с помощью критерия Фишера Р н 
по условию целесообразности: Р н >Г( Р.к ь к 2 ). 


Поскольку для химико-технологических расче- 
тов доверительная вероятность Д=0,95 [9]; и для 
данной модели числа степеней свободы к х = 3, к 2 = 3, 
то ДДЛЛ)=9,3 [10]. 

Параметры разработанной модели представле- 
ны в табл. 3 и свидетельствуют о ее высокой эффек- 
тивности. 


Таблица 3. Параметры разработанной модели 


Параметр 

Флегмовое число | 

1 

2 

3 

4 

У, 

45,97 

43,82 

42,95 

42,49 

Л а 

20,67 

0,00 

3,43 

8,00 

5 У 

2195,82 

2260,64 

2294,44 

2314,27 

Ру 

106,22 

5023634,96 

668,50 

289,28 


у - среднеарифметическое всех измерений; 5% - дисперсия аде- 
кватности; $ - дисперсия относительно среднеарифметического 


Традиционно оценка пригодности контактного 
устройства для проведения процесса определяется 
по величине межтарельчатого уноса е, оптималь- 
ная величина которого е ор , находится в интервале 
[0,2. ..0, 32] кг жидкости/ кг пара. 

На основе текущих составов потоков и их фи- 
зических свойств получены его значения по высоте 
колонны (рисунок). 
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Рисунок. Распределение межтарельчатого уноса по высоте 
ректификационной колонны при различных вели- 
чинах орошения (г - номер ступени разделения) 


Очевидно, что в колонне только при минималь- 
ном флегмовой числе значения межтарельчатого 
уноса приближаются к оптимальным, и лишь в ук- 
репляющей части. В целом процесс в ректифика- 
ционной колонне проводится при высоких величи- 
нах межтарельчатого уноса жидкости, что свиде- 
тельствует о необходимости модернизации дей- 
ствующих контактных устройств. 
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Таким образом, выполненная работа позволяет 
моделировать процесс промышленной ректифика- 
ции смеси Ф-М-В, осложненный химическими 
взаимодействиями между компонентами, адекват- 
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или иных типов контактных устройств, оптимизи- 
ровать технологический режим работы колонны. 
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